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ABSTRACT 

Observations and modelling suggest that the fluctuations in magnetised plasmas exhibit scale- 
dependent anisotropy, with more energy in the fluctuations perpendicular to the mean magnetic field 
than in the parallel fluctuations and the anisotropy increasing at smaller scales. The scale-dependence 
of the anisotropy has not been studied in full-orbit simulations of particle transport in turbulent plas- 
mas so far. In this paper we construct a model of critically balanced turbulence, as suggested by 
iGoldreich fc Sridhaii ()1995l ). and calculate energetic particle spatial diffusion coefficients using full- 
orbit simulations. The model uses an enveloped turbulence approach, where each 2-dimensional wave 
mode with wavenumber k± is packed into envelopes of length L following the critical balance condi- 
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tion, L (xk^ ' , with the wave mode parameters changing between envelopes. Using full-orbit particle 
simulations, we find that both the parallel and perpendicular diffusion coefficients increase by a factor 
2, compared to previous models with scale-independent anisotropy. 
Subject headings: cosmic rays - diffusion - turbulence 



1. INTRODUCTION 

The origin of Solar Energetic Particles (SEPs) is one 
of the unsolved problems in heliospheric physics. Both 
flares and coronal mass ejections are capable of ac- 
celerating particles, and studies analysing SEP events 
give differing conclusions on the main particle accelera- 
tors (see, e.^ . ICane et al.ll201CII : IGopalswamv et aLll2012l : 
lAschwa nden 2012). The interpretation of SEP events 
is complicated due to the particle propagation effects 
caused by the turbulent interplanetary magnetic field. 
SEP events are obse rved at a large range of latitudes 
and longitudes (e^ ., iDalla et all 120031 iLiu et al.l 120111 : 
iDresing et al.112012 1. suggesting possible cross-field trans- 
port of the SEPs. Thus, in order to understand the SEP 
acceleration in the solar eruptions, we must first under- 
stand the nature of SEP tranport in the heliosphere. 

Cosmic ray research is typically based on two the- 
oretical appro aches: the diffusion-convection equation 

( 

( 
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and the quasilinear approach (QLT) 
The original description of the cross-field 
transport of a charged particl e as propagat ion in random- 
walking magnetic field lines pokipiilll966[ ) has been im- 
proved to take into account the interpla y between parallel 
and perpendi cular propagation effects (jMatthaeus et al.l 
l200alShalchill2010f l. The spread of the particles observed 
in the SEP events, however, remains difficult to explain, 
as they require larger ratios of the perpendicular-to- 
parallel diffusion coefficients tha n what can be expec ted 
from the theoretical approaches (jDresing et al.|[201^ . 

Recently, particle transport research has benefited 
from numerical simulations utilising particle full-orbit 
simulations, which have the advantage of needing 
no a priori assumption s for the particle prop a gation . 
iBeresnvak et al.l (|2011[ ) and iWisniewski et al.l (|2012[ ) 
have recently studied particle propagation using MHD- 
simulated turbulence. However, as MHD simulations are 
limited to small range of scales, most full-orbit parti- 
cle simulations describe the fiuctuating fields as a super- 
position of Fourier modes on a constant magnetic field 



e.g.lGiacalone fc Jokipiilll999l:lOinll200l[Qin et al.| [2002l: 
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Observations and modelling suggest that the fiuc- 
tuations in magnetised plasmas are anisotropic, with 
more energy in the fluctuations perpendicular to the 
mean magnetic field than in the parallel fiuc tuations 
(e.g. IShebalin et all [19831: IBieber et al.l 119961) . The 
anisotropy is also scale-dependent, with structures at 
smaller scales more anisotropic than the larger scales. 
The scale-depende nce o f the a nisotropy was predicted 
by Goldreich & Sr idha^ ([19951 hereafter GS95), who 
noted that as the turbulence amplitudes increase at small 
scales, they introduce variation along the mean field di- 
rection. This is due to the fact that when interacting, the 
wave packets propagate through each other, along the 
mean magnetic field. When the changes in the interact- 
ing wave packets reach nonlinear amplitudes, both of the 
interacting waves change, and for this reason also the out- 
come of the interaction changes. This takes place on non- 
linear timescales tmli when the wave packet's front has 
propagated distance oc Va t^l, where Va is the Alfven 
velocity. Using the nonlinear timescale from Kolmogorov 

2/3 

scaling, GS95 obtained a critical balance, A:|| oc k_^ 
between the parallel and perpendicular wavenumbers. 
The critical balance between structure scales has been 
found in MHD turbulence simulations (e.g., Clio et al.l 
'2002*) , and observed in the solar w ind turbulence (e.g., 
■Horburvet all 120081 : iPodestal [20091 ). In order to under- 
stand the propagation of particles in heliospheric plas- 
mas, the scale-dependence of the turbulence anisotropy 
should thus be taken into account. 

In this work, we study the particle propagation in 
the critically balanced turbulence of GS95. To model 
the turbulenc e, we use the enveloping approach of 
iLaitinen et al.l ()2012D . The enveloping approach im- 
proves the earlier models that use infinite, linear plane 
waves by breaking the phase coherence of the waves in 
the direction of the mean magnetic field. Such a loss of 
coherence is expected in the plasma in the heliosphere. 
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as the waves evolv e in a nonline ar manner (se e , e.g. , 
iTu fc Marsc"hlll995l for a review). iLaitinen et all (|2012[ ) 
studied the effect of structures on the tranport and used 
a fixed envelope length of the order of turbulence cor- 
relation length. In the present paper, we study the 
particle propagation in turbulence with scale-dependent 
anisotropy, and the enveloping is done separately for each 
wave mode, with the envelope length following the crit- 
ical balance scaling, L oc fc^ , where L is the enve- 
lope length and k± the wavenumber of the enveloped 
Fourier mode. The Fourier modes are chosen to be 2D 
modes, with their wave vector perpendicular to the mean 
magnetic field. Thus, the variation of the magnetic field 
along the mean field is purely due to the enveloping of 
the modes. We use particle full-orbit simulations in the 
modelled turbulence, and calculate the diff'usion coeffi- 
cients for energetic particles, comparing them to results 
given in previous work. In Section 2 we describe the crit- 
ically balanced turbulence model and the method used 
for obtaining the energetic particle diffusion coefficients. 
In Section 3, we study the scale-dependence of the tur- 
bulence in the developed model, and report changes in 
the energetic particle diffusion coefficients, compared to 
the traditional composite model. We discuss the results 
in Section 4, and draw our conclusions in Section 5. 

2. MODEL 

2.1. Turbulence Model 

In our model, the magnetic field consists of a uniform 
and constant background field, Bq, overlaid with a turbu- 
lent field, (5B(x, y, z), with the total magnetic field given 
as 

B = Boe^ +SB{x,y,z), 

For the background field we use Bq — 5 nT, consistent 
with the magnetic field magnitude at 1 AU. 

We model the scale-dependent turbulence by us- 
ing the turbu l ence envelope approach, introduced by 
ILaitinen et al.l ()2012l ) . In the approach, the turbulence is 
enveloped into packets along the mean field direction (see 
Fig.[T]). From an envelope to the next, the parameters of 
the wave modes (random phase, wave vector direction) 
change. Thus, a particle propagating from an envelope 
to another interacts with a coherent Fourier mode only 
within one envelope, and the linear coherence is broken 
when it enters a different envelope. 

The scale-dependence of the turbulence in our model 
is achieved through the selection of the envelo pe lengths, 
which in the study of ILaitinen et al.l ()2012D were con- 
stant. Here the enveloping is done separately for each 
wave mode, with the wave scale Aj_„ = 2Tr/k±„, where 
k±n is the wavenumber of the 2D mode. The enve- 
lope length, Ln, follows the critical balance scaling, 

Ln = Lij Xj^^ , where Lc is the scale where parallel and 
perpendicular scales are in balance. As a result, for wave 
scales Aj.„ < Lc the perpendicular wave scales decrease 
faster than the parallel envelope scales, resulting in scale- 
dependent anisotropy. This process is depicted in Fig. [H 
where, we show the wave scale (thick horizontal line) 
compared to the envelope length scales, for four different 
wave scales. 
The amplitude of mode n for envelope i is modulated 
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Fig. 1. — The green and white packets represent the envelopes 
along the mean field direction, for four different perpendicular 
scales, shown by the thick horizontal lines. The steepness param- 
eter is cr = 10 in this figure. 
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where Zi marks the beginning of the envelope and z is 
the distance along the direction of the constant magnetic 
field. The profile of the envelope is given by function 



F{z} 
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z<-l 
-1 < z < 1 

z> 1 



(1) 



This approximates the often-used tanh profile, which en- 
ables analysing differentiable profiles that can have steep 
edges, approaching a step-function shape. The polyno- 
mial description was chosen for numerical eflficiency. The 
steepness of the envelope profile is given by parameter 
Sn = Ln/a, and cr is is a dimensionless parameter relat- 
ing the envelope length to the rise profile length. 

In the modelling, we consider that a wave mode with 
a given wave vector direction and random phase gives its 
energy to another mode with the same modulus of fc„, 
without losses, thus keeping the sum of the amplitudes 
of the modes constant, J2i^n,i{z) — 1. This is achieved 
by selecting z^+i = Zi + Ln- 

The steepness of the envelope is related to the trans- 
fer rate of energy from a mode to another, due to the 
critical balancing of the turbulence. The rate of change 
depends on the relative wavenumber directions of the in- 
teracting modes (e.g.Luo & Melrose 2006), which results 
in some interactions transferring the energy faster than 
others. In our model, however, we describe phenomeno- 
logically only the result of the critically balanced 3- wave 
interactions, thus accurate description of the rise pro- 
file is not possible on this level. Instead, we consider in 
this study two different profiles, a gradual change with 
a — 1 and a rapid change, with cr = 10. As shown in 
Fig. [21 for cr = 10 the energy is transferred from one 
mode to another in an intermittent manner, whereas for 
cr = 1 the interaction is more gradual, between multiple 
modes. The distance along the mean field direction over 
which one mode dominates over the others, remains the 
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Fig. 2. — The green and white packets represent the envelopes 
along the mean field direction, for cr = 10 (top panel) and a = 1 
(bottom panel). The envelope length L„^i = 1 for these profiles. 

same, the envelope length, ij„, for both values of a. 

The turbulence field is given as a sum of Af„ envelopes 
and N Fourier modes, as 



dB{x,y,z) 



N M„ 
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In this study the Fourier modes in the envelopes are 
2D modes, with the wave vector k_L„ perpendicular to 
the mean magnetic field. Thus the variation on parallel 
scales in the model are due to the enveloping only. The 
polarisation vector ^„^i lies in the xj/-plane and is per- 
pendicular to both the mean magnetic field and k„, in 
order to satisfy V • B = 0. Thus 
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where the coordin ate system r^ is obtained from r with 
a rotation matrix (jGiacalone fc Jokipiilll999l ). The wave 
propagation directions in the xy-plane, as well as the 
random phases /3„^i, are chosen from a uniform random 
distribution. 

The fluctuation amplitude A{kn) is given by a power 
law spectrum 
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where Bf is the variance of the magnetic fleld, Lc the 
spectrum's turnover scale, 7 the spectral index, and AT/J^ 
specifles the volume element in fc-space that the dis- 
crete mode k±n represents. For the turnover scale we 
use Lc = 2.15 r0, where r© is the solar radius. For 
2D turbulence, the spectral index is 8/3 and the factor 
AVn = 27rfc„Afc„. In the simulations we use logarithmi- 
cally spaced wavemodes with the wavenumber running 
from 27r/lAU to 27r/10-'^AU. 

To compare the turbulent fields and their effects on 
SEP propagation, we scale the fluctuation amplitudes so 
that the average energy density in the fluctuating fleld is 
independent of the enveloping parameters. The scaling 
factor is obtained through numerical integration. 

2.2. Energetic Particle Simulations 



We study particle propagation in the modelled turbu- 
lent magnetic flelds by integrating the fully relativistic 
equation of m otion of en ergetic protons using the simu- 
lation code bv lDalla fc B rowning (2005). The code uses 
the Bulirsh-Stoer method (Press ct al.l ll993l ). with adap- 
tive timestepping to control the accuracy by limiting the 
error between consecutive steps to a given tolerance. We 
simulate 2048 particles in 10 fleld realisations, with dif- 
ferent wave mode random phases and wave vector di- 
rections in each realisation, thus giving a total of 20480 
particles. Each realisation has N=128 Fourier modes, 
whereas the number of envelopes depends on the energy 
of the particle. The particles are simulated for ~ 100 
parallel diffusion times. The particle diffusion coefficient 
is obtained with 
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(e.g. iGiacalone fc JokipiJlll999D . with kh equal to k^ 



and 
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whereas kj_ is obtained as the mean of Kx. 
further details of the simu lation scheme, and v erifica - 
tio n against the results of IGiacalone fc Jokipnl (J1999D 
see lLaitinen et al.l ()2012[ ). 



3. RESULTS 



3.1. 



Scale- dependence of the turbulence 

The goal of this work is to study particle motion in 
a turbulent field that has a scale-dependent anisotropy 
that corresponds to the critically balanced scaling sug- 
gested by GS95. As the fi rst step, we compare the co m- 
posite turbulence model of lGiacalone fc Jokipiil (|1999D to 
the one presented in this study. We do this by calculating 
the two-point correlation function. 



CxA^)^{SBx(r)6Bx{r + s)) 



(4) 



where s is the scale vector the correlation is calculated 
for, and the angle brackets denote ensemble average. We 
calculate the correlation function as a function of s^ 
and S||, that is, across and along the mean field direc- 
tion. The correlation function is calculated from the dis- 
sipation scales to the turnover scale Lc of the spectrum, 
and presented in Fig. [3l where the contours are selected 
to represent correlations at equidistant scales, s, rather 
than values of the correlation function, Cxx ■ The correla- 
tion contours for the GS95 model are ellipses (appearing 
rectangular in the log- log presentation of Fig. [3]), with 
aspect ratio increasing at smaller scales. A similar trend 
cannot be seen in the composite model. 

To quantify the scale-dependence of the anisotropy, we 
calculate the axis ratio of the contours. This is done by 
finding values of the parallel and perpendicular scales, 
S|| and s±, that satisfy Cxxis±,0) — Cxx{0,s\\). We 
present the axis ratio sh/s± as a function of sj_, for the 
turbulence models, in Fig. SI The dash-dotted straight 
fine represents the GS95 scaling, the solid and the dot- 
ted curves (green and red in the online version) repre- 
sent the phenomenological model of this work, with two 
different values of the steepness parameter ct, and the 
dashed curve (blue in the online version ) repre sents the 
composite model of lGiacalone fc Jokipiil (|1999D . As can 
be seen, the axis ratio in the model presented in this 
work is clearly scale-dependent, and follows the trend of 
the GS95 scaling. There is a significant dependence of 
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Fig. 3. — Two-point correlation functions for the composite 
model (top panel), and for the GS95 turbulence model (bottom 
panel). 



the level of anisotropy on the steepness parameter, with 
the steep profile for the energy change between the wave- 
modes, represented by cr = 10, producing values closer 
to isotropy when the axis scales approach the isotropy 
scale Lc- 

At large scales, the phenomenological model deviates 
from the GS95 model. This is due to the different form 
of the turbulence spectrum used in the derivation of 
the critical balance by GS95, and the one used in this 
study. The GS95 scaling is obtained by using a con- 

5/3 

tinuous power law spectrum, k^ , representing the Kol- 
mogorov turbulence inertial range spectrum, while in the 
heliosphere the inertial range is limited in extent to only 
a few orders of magnitude, as the spectrum flattens at 
large scales into energy-containing range. This results 
in flattening of the correlation function at scales of or- 
der Lc and larger. Thus, when Aj^ < L|| < Lc, with 



the parallel scale L|| reaches the flatten- 
ing of the correlation function before the perpendicular 
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Fig. 4. — The scale-dependence of the two models with the axis 
ratio su/sj_. In the grey area, LJ Sj^ exceeds the turnover scale 

scale A_L, and the ratio Lu/X± remains anisotropic. This 
range is shown by gray shading in Fig. 21 As the spec- 
trum deflned by Eq. (jS]) rolls gradually to the inertial 
range, it is expected to see the deviation from the GS95 
form to start already at somewhat smaller scales. Solar 
wind turbulence observations show a somewhat smaller, 
but still anisotropic axis ratio of 1.5 for scales of ~ tq 
for slow solar wind ()Dasso et al.ll2005l ). 

Thus, we conclude that our model with cr = 10 dis- 
plays the required scale-dependence of anisotropy, with 
the large scale anisotr opy being co n sisten t with the solar 
wind observations of iDasso et al.l ()2005D . The gradual 
energy transfer model, a = 1, produ c es hig her level of 
anisotropy, exceeding the IDasso et al.l ()2005[ ) result. 

3.2. Turbulence spectrum 

In the previous section, we described the scale- 
dependent behaviour of the two-point correlation func- 
tion for different turbulence models. When studying tur- 
bulence, it is often more common to present a power 
spectrum, which quantifies the energy deposited in the 
turbulence at different scales. We present such spectra in 
Fig. [5l The spectra are obtained from one-dimensional 
two-point correlation functions, with parallel spectrum 
defined as 



L^xxyf^W 
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and the Pxx{k±) similarly from Cxx{s±,(^)- The one- 
dimensional spectrum can be compared to the three- 
dimensional spectral form given by GS95 as 
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where C is a dimensionless constant, L the isotropic ex- 
citation scale, and g{x) a function that vanishes at large 
X. From this form the one-dimensiona l spectrum c an be 
obtained using the Taylor hypothesis (|Tavlodll938| ). 
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Fig. 5. — Parallel (black curve) and perpendicular (grey curve, 
green in the online version) power spectrum of the GS95 turbulence 
model with a = 1 (top panel) and tr = 10 (bottom panel). The 
dash-dotted and dashed curves depict the theoretical spectra given 
by Eqs. [5] and [B] respectively. The vertical bars at the top of the 
panels depict the inverse of the larmor radii of protons of energy 
100, 10, 1 and 0.1 MeV, respectively, from left to right. 



ft is straigthforward to show that the one-dimensional 
spectrum along the field line is proportional to fc~^ and 
the perpendicular spectrum to k^^'^ . Furthermore, we 
find that when using the simpler form suggested by 
iCho et all ()2002l ). g{x) = exp(— x), the spectra are given 
in form 
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We plot these spectra in Fig. [SI using Lc for the isotropic 
excitation scale L, and fitting the perpendicular spec- 
trum, Eq. ^ , to the perpendicular spectrum of the sim- 
ulated turbulence. It should be noted that the parallel 
spectrum represents a continuum instead of being com- 
posed of discrete modes, as it results from the envelop- 



ing rather than begin generated with a sum of parallel 
Fourier modes. 

As seen in the figures, the envelope shape affects the 
parallel spectrum shape. For the gradual envelope pro- 
file, with a = 1 (Fig.[5l top panel), the parallel spectrum 
has spectral index of -2.5, steeper and an order of mag- 
nitude below the spectrum given by Eq. ([S]). For the 
steep profile, with a = 10 (Fig. [Sj bottom panel), the 
parallel spectrum has a spectral index of —2, consistent 
with Eq. ([SJ, with the relative power in the parallel and 
perpendicular spectra consistent with Eqs. ^ and ^. 
The perpendicular spectrum is not affected by the en- 
veloping, with the spectral index remaining at the Kol- 
mogorov value of -5/3, the input of the model for both 
values of a. The spectral index of -2 for the parallel 
direction and -5/3 for the perpendicular direction were 
recent l y observ ed in the solar wind by iHorburv et al.l 
([200l IPodestal (2009) . The ratio P±/P\\ in our model is 
2 at fc = 10 r^^ and 5.7 at 100 rZ^, which is similar to 



the results of IHorburv et al.l ()2008l ) and IPodestal (|2009D . 
who find values of 1.5-5 for the range k = 10 — 1000 Tq^. 

3.3. Energetic Particle Propagation 

In order to understand how the scale-dependence of 
turbulence anisotropy affects the energetic particle prop- 
agation, we have simulated energetic protons at energies 
0.1, 1, 10 and 100 MeV, with the particles Larmor ra- 
dius scale, fci = r£ shown as vertical bars in Fig. \E\ 
To study the collective behaviour of the particles, we 
have calculated the particle diffusion coefhcients along 
and across the mean magnetic field, as described in Scc- 
tion l2.2[ shown in Fig.[6l with statistical error limits. For 
comparison, we have also followed particles in the com- 
posite model turbulence of Giacalone & Jokipii (1999). 
The parallel diffusion coefficients for the GS95 model and 
the composite model are presented in the top panel of 
Fig. ini In the figure, the solid curve represents the com- 
posite model, whereas the dashed and dash-dotted curves 
(green and red in the online version) represent the GS95 
model, for two values of the steepness parameter. As can 
be seen, the parallel diffusion coefficient is larger for the 
GS95 model, and depends strongly on the steepness pa- 
rameter of the scale-dependent enveloping. This can be 
understood on the basis of lower energy in parallel fluctu- 
ations to scatter the particles in the a — 1 case compared 
to the cr = 10 case. The dependence of the parallel dif- 
fusion coefficient on particle energy in the GS95 model 
also differs from the composite model, with the diffusion 
coefficient decreasing slower at decreasing energies. This 
can be understood as the effect of the steeper spectrum, 
which results in smaller relative energy in the small-scale 
fluctuations that scatter the lower-energy particles. 

We present also theoretical estimations of the parallel 
diffusion coefficient in F ig. El The dotted line repre sents 
the QLT K||, as given bv lGiacalone fc Jokipiil ()1999f ) pre- 
sented in their work for comparison with the compos- 
ite and isotropic turbulence models. Their estimation 
is based on an equation that is valid only for parallel 
Alfven waves (the slab waves). It should be noted that 
the in a composite turbulence the 2D component does 
not have a significant cont ribution to parallel scatter- 
ing (e.g.. lBieber et anil994D . If we assume that only the 
slab component contributes to the parallel scattering, the 
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coefRcient can be estimated by 
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Fig. 6. — Parallel (top panel) and perpendicular (bottom panel) 
energetic particle diffusion coefficients in composite turbulence 
(solid curve), and in the GS95 turbulence (dashed and dash-dotted 
curves, green and red in the online version), as a function of proton 
energy, for turbulence amplitude Bi/Bq = 1. In the top panel, the 
dashed curve depicts the QLT diffusion coefficient for the composite 
model, as given bv fGiacalone fc Jokipii ( 1999i ), and the tripledot- 
dash line Eq. [7]for the a = 10 case. 

80%:20% mix of 2D and slab modes implies a factor 5 rise 
of the dotted curve in the upper panel of Fig. |6l which 
fits well to the simulated composite model results. 

Comparison of the theoretical parallel diffusion coef- 
ficient with our results for the GS95 turbulence is not 
straightforward, as the QLT for magnetostatic turbu- 
lence results in infinite diffusion coefficient for parallel 
spectrum ex ^7 and steeper. This resonance gap prob- 
lem, due to lack of scattering at low particle pitch angle 
cosines, fi — v\\/v, has been studied in the scientific liter- 
ature, and seve ral mechanisms for its closure have been 
suggested (e.g., 'Schlickciser & Achatz"1993|; iBieber et al.l 
11994 iNgfc Reanies 1995; Vainio 2000). As shown by, 
e.g.. lVainiol (|2000f )7if scattering at low pitch angles after 
the closure of the resonance gap is weak, the diffusion 
coefRcient is dominated by this low rate of scattering, 
and thus by the closure mechanism. If, however, the res- 
onance gap is efficiently closed, the pitch angle diffusion 
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where r^ is the particle's Larmor radius, and h{fj,g) is 
a dimensionless function that depends on the spectral 
index of the parallel turbulence and the pitch angle fig 
below which the resonance gap closure mechanisms over- 
come the QLT result. This form is valid only for a power 
law spectrum of the turbulence, and thus the higher en- 
ergies in this study are expected to have larger diffusion 
coefficients, compared to the estimation. In addition, it 
should be noted that the above considerations are valid 
for slab turbulence. As our model contains a wealth of 
oblique wave modes, implied by the correlation function, 
the QLT would give a wealth of harmonic resonances in 
addition to the resonance at inverse Larmor radius. 

We present the theoretical parallel diffusion coefficient, 
as given by Eq. ([7|), in Fig. |6l with a tripledot-dash line, 
for our model with steepness parameter a = 10. The 
power spectrum in Eq. ([7]) is obtained as a power law 
fit to the parallel spectrum (Fig. [5]). The exact scaling, 
given by h(fig), is not known, thus we can only com- 
pare the trends of the diffusion coefficients. As can be 
seen, the QLT trend tracks the GS95 model only at low 
energies. This suggests that the resonance gap effects 
must be taken into account when estimating the mean 
free path at steep spectra. 

The perpendicular diffusion coefficients are shown in 
the bottom panel of Fig. [6l The perpendicular diffusion 
is stronger in the GS95 turbulence, by a factor of 2, as 
compared to the composite model. The dependence on 
the enveloping shape is small except for the highest en- 
ergies used in the study. Thus, the perpendicular and 
parallel diffusion coefficients behave differently with re- 
spect to the exact description of the critically balanced 
turbulence. 

4. DISCUSSION AND CONCLUSIONS 

In this work, we study how energetic particles prop- 
agate in critically balanced turbulence. We build tur- 
bulence with scale-dependent anisotropy resulting from 
critical balance by using envelopes of length that fol- 
lows the required scaling between the perpendicular 
wavenumber and the parallel variation scale. Analy- 
sis of the 2-point correlation function and the spectra 
shows that for the envelope steepness parameter cr = 10 
our model is consistent with the GS95 critical balance 
scaling and the observations of solar wind turbulenc e 
(|Dasso et al.l [20051 : iHoTburv et al.l[200l iPodestal [2009( ) . 
The gradual shape, cr = 1 (see Fig. [2]) results in more 
pronounced anisotropy. The case of cr = 10, however, 
can be considered more realistic, as the nonlinear in- 
teractions be tween the waves dep end on their relative 
phases (e.g., iLuo fc Melrosel[2006[ ). resulting in an un- 
even evolution of the wave packets. The structure of 
rapid change after an extended period of unchanged tur- 
bulence can also be viewed as intermittent turbulence. 
We note, however, that in our model the probability dis- 
tribution of magnetic field increments, ABx{s±, sn) — 
Bxiri_ + s_L-''|| + s\\) ^ Bx(r±,rn) does not significantly 
deviate from Gaussian , unlike fo r the observed and mod- 
elled turbulence (e.g.. , Greco et al.,,2009.) . 
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The scale-dependence of the turbulence results in a 
notable difference in the power spectra along and across 
the mean magnetic field, as shown by Fig. [5] Thus, it 
is expected that the energe tic particle transport differ s 
from the composite model of lGiacalone fc Jokipiil (J1999D . 
where the spectral indices are both equal to 5/3. Our 
simulations show that the effect for the parallel diffusion 
coefficient is significant, with a 1.5-2.5 -fold increase, de- 
pending on the particle energy, for cr := 10 (top panel of 
Fig. [6]). An increase of the diffusion coefficient with a 
steepening parallel spectrum is consistent with the stan - 
dard quasilinear theory (e.g.. lDung fc SchlickeiseririQQOl ). 
Our analysis shows, however, that the simple estimate 
based on QLT, for the parallel wave modes (Eq. dZ])) fails 
to describe the parallel diffusion coefficient at higher en- 
ergies. 

The perpendicular propagation is also affected by the 
scale-dependence of the turbulence, with the perpendic- 
ular diffusion coefficient larger by a factor 2 compared 
to the one obtained with composite m odel. This can be 
under stood on the basis of the work of iMatthaeus et all 
(|2003l). In their formulation, the perpendicular diffusion 
coefficient depends on the parallel diffusion coefficient. 
Thus, a different form for the parallel component spec- 
trum would influ ence the cross-field tr ansport of the en- 
ergetic particles. iShalchi et al.l (120101 ) used the further- 
developed model bv lShalchil ()2010[ ) to calculate the per- 
pendicular diffusion coefficient in a GS95-type spectrum. 
Their kj_/k|| decreases '~ 5-fold compared to the com- 
posite model, and deviates from our result, where the 
ratio does not change for a — 10. However, in their work 
they did not estimate the parallel diffusion coefficient 



consistently using the GS95-spectrum , but used stan- 
dard quasilinear theory pokipiilll966f) with Kolmogorov 
spectrum, and a model based on Galactic cosmic ray ob- 
servations. Thus, their results are not comparable to 
k^/k|| in our model. The insensitivity of k±^ to the par- 
allel diffusion coefficient in their model is in line with 
our results of insensitivity of k_l to the changes in the 
turbulence caused by the change of the envelope shape. 
This inse nsitivity can be implied also from the theoretical 
work bv IMatthaeus et all (|2003[ ) (their Eq. (7)), where 
the rapidly decreasing power spectrum may diminish the 
link between the parallel and perpendicular diffusion. 

Determining the correct scattering parameters in the 
interplanetary space is very important for the analysis 
of SEP events. Studies have shown that the scatter- 
ing along the interplanetary magnetic field can signifi- 
cantly affect the accuracy of the SEP event onset analy- 
sis a t 1 AU (iLintunen fc Vainiol [200l iSaiz et al.l [20051: 
Laiti nen et al.l I2O1O0 . Recentlv. iGiacalone fc Jokipiil 
(|201^ showed that varying the diffusion coefficients and 
the ratio of the perpendicular and parallel diffusion co- 
efficient in the simulations can result in a wide vari- 
ety of intensity profiles at 1 AU for an impulsive SEP 
event. Therefore the two- fold increase in the parallel and 
perpendicular transport coefficients reported in this pa- 
per will significantly influence the characteristics of SEP 
events. 
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